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DUST ECHOS FROM GAMMA RAY BURSTS 
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ABSTRACT 

The deviation from the power-law decline of the optical flux observed in GRB 970228 and GRB 980326 
has been used recently to argue in favor of the connection between GRBs and supernovae. We consider 
an alternative explanation for this phenomenon, based on the scattering of a prompt optical burst by 
O.IMq of dust located beyond its sublimation radius 0.1 — 1 pc from the burst. In both cases, the optical 
energy observed at the time of the first detection of the afterglow suffices to produce an echo after 
~ 20 — 30 d, as observed. Prompt optical monitoring of future bursts and multiband photometry of the 
afterglows will enable quantitative tests of simple models of dust reprocessing and a prediction of the 
source redshift. 

Subject headings: gamma rays: bursts - ISM: dust, extinction 



1. INTRODUCTION 

The relationship between Gamma Ray Bursts (GRBs) 
and supernovae has become increasingly interesting over 
the past year. Though exploding massive stars have long 
been consider ed as possible progenitors of GRBs (e.g. 



Woosley 1993 ), no evidence existed to support these the- 



ories until observations of the afterglow of GRB 980425 
suggested an assoc iation of the burst with an unusual 



super nova 1998bw ( Galama et al. 1998 , Kulkarni et al 



1998 ). Later reanalysis of the optica l afterglow lightcurves 
of two other burst s, GRB 970228 flFruchtcr et al. 1999[ ) 



and GRB 980326 flBloom ct al. 1999j ) showed a deviation 
from the power-law decline expected if the emission is due 
to synchrotron radiation from electrons accelerated by the 
blast wave. In both cases a significant excess emission was 
observed around ~ 30 days after the gamma- ray burst, 
with simultaneous reddening of the spec trum. |Bloom ct 



al. (1999)| , |Galama et al. [1999J and |Reichart (1999) 



attribute this excess to the emission from an underlying 
supernova event. 

The relationship of GRBs to SN explosions is a question 
of great importance, since it provides a powerful clue to 
the fundamental nature of these objects. However, the evi- 
dence presented so far is circumstantial - the association of 
GRB 980425 with SN 1998bw is unproven and the excess 
emission seen from GRB 970228 and GRB 980326 is based 
upon relatively few actual measurements - and possible al- 
ternative explanations need to be seriously considered, if 
only t o strengthen the case for th e SN explanation. In this 
spirit, Waxman fc Draine (1999) suggested that the red ex- 
cess emission observed in GRB 970228 and GRB 980326 
is due to dust in the vicinity of the burst progenitor ab 
sorbing and then re-radiating the optical/ UV flash ob 



serve d shortly after the recent GRB 990123 ( Akcrlof et al 



1999 ) and generally attributed to t he reverse shock which 
propagates into th e fireball ejecta flMeszaros, Rees fc Pa 



pathanassiou 1994 , Meszaros fc Rees 1997 , Panaitescu <L' 



Mcszaros 1998, Sari fc Piran 1999). However, the Wax- 



man & Draine scenario has two shortcomings. Firstly, the 
equilibrium temperature of dust is limited to ~ 2300 K 
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and so the emission should peak at ~ 2 (1 + z) /zm (where 
z is the GRB redshift), although a small amount of higher 
temperature emission may be produced by the dust as it is 
subliming. Secondly, the optical flash is so powerful that 
the sublimation radius lies beyond ~ 10 pc from the GRB. 
Thus, in this picture it is rather difficult to reproduce the 
observed flux in the 0.4 — 0.8 (im band with a time delay 
of order a few weeks. 

In this letter we propose an alternative explanation, 
which relies on the scattering of the direct optical transient 
emitted in the first day by dust as the primary source of 
excess optical radiation. The fundamental point is that in 
the two observed cases, assuming isotropic emission, the 
fluence of the observed transient exceeds that of the re- 
ported excess and the unobserved transient is even larger 
if we extrapolate to earlier times. A fraction of this emis- 
sion scattered from a radius where dust can outlive the 
optical transient should therefore produce a delayed echo. 
As dust absorbs selectively as well as scatters, the echo is 
likely to be significantly redder than the original optical 
transient, as reported. 

In the next section, we describe our model for the dust 
scattering properties and then present the results in the 
context of the observed GRBs in §||. Implications for fu- 
ture tests of our scenario are discussed in §|]. We assume 
h = 0.6, r^M = 0.3, and ft a — 0.7 so that the angular 
diameter distance of the GRB is Da = 1.5 — 2 Gpc for 
0.5^z^3. 

2. DUST ECHOS 
2.1. Sublimation Radius 



Waxman & Drainc (1999) estimate that dust grains in 
the path of the optical/UV flash will be effectively sub- 
limed out to a distance 



■Rsub ~ 1 (Qab s £47/a-i) 1/2 pc, 



(2-1) 



where Q a bs ~ 1 is the absorption efficiency factor for op- 
tical/UV photons, L4J = J dvL u /lQ i7 erg s -1 is the un- 
beamed luminosity of the optical transient (OT) in the 



2 



1 — 7.5 eV energy band, a_i ~ 1 is the dust grain size in 
units of 0.1 fjxa. Beyond R su b, only the most refractory 
grains, like silicates, can survive. Note that the thermal 
time for a typical dust particle is of order 1CP 4 — 10~ 2 s, 
much shorter that the duration of the optical transient, 
so that we can treat grains as being in thermal equilib- 
rium with the incident radiation, (which has a pressure 
-Psub ~ 0.03 dyne cm~ 2 at i? S ub)- 

The extinct ion properties of silicat e dust particles were 
computed by Draine & Lee (1984) (see their Fig. 10), 



for a power-law distribution of part icle sizes proposed by 
Mathis, Rumpl fc Nordsieck (1977)| to explain interstellar 
starlight extinction. Based on their results, we take the 
ratio of the scattering and absorption efficiency factors to 
be of order Q sc /Qabs — 4, and the average scattering an- 
gle to be (cos0) = (/u) ~ 0.5 for observed wavelengths 
0.2- 1 (l + z)/im. 



~ Ai OT . We expect that, in practice, it will be quickly 
destroyed by the effects of secondary cosmic ray electrons 
created through electron scattering of the GRB so that the 
observed optical afterglow need not necessarily be subject 
to the same extinction as the echo. 

2.3. Optical Scattering 

The optical echo flux density , observed at frequency 
is oh = + z)- 1 is 



^f ah { v oh,t a b) 



L v {y,t) cAt OT dV sc 
2(1 + z) 3 D\ R dTi 

0<t oh < 2R suh (l + z)/c, 



(i/,/*); (2-3) 



where dV sc (v, /x)/dfi, is the probability of escape along 
the direction defined by angle 6 — cos -1 [i for a photon of 
frequency v. 
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Fig. 1. — A schematic diagram of the GRB environment. The 
long-dashed line represents the position of the expanding optical/UV 
photon front at time t in the frame of the GRB. The shaded area 
shows the region where the dust is not sublimated instantaneously. 
The heavily shaded area shows the region of the shell from which the 
scattered radiation is observed while the lightly shaded area repre- 
sents the regions where the dust is not scattering any more. 
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Fig. 2. — Solid and long-dashed lines show the escape probability 
for photons scattered by a dust slab for different values of 9 and r 
(as marked in the figure). For comparison, the dotted line represents 
an escape probability that increases the spectral index of the echo 
relative to the OT by 2. The thin long-dashed line shows the results 
computed using a different dcr/d£l oc 1 + 2fi + fi 2 which gives similar 
results although it is less peaked at fi ~ 1. The short-dashed line 
shows the fraction of photons at each wavelength, /ns(A) which pass 
through the slab unscattered, approximating the escape probability 
at 6 = 0. 



2.2. Source Geometry 

Fig. [I] shows a schematic picture of the GRB environ- 
ment observed at time t after the detection of 7-rays. The 
incident optical transient emission is supposed to be lim- 
ited to an interval At OT = Ai° b T /(l + z) « t after the 
GRB, and scattered by dust beyond i? su b- We specialize 
immediately to the case when the dust is associated with 
an outflowing spherical wind, and the OT is isotropic. (It 
is straightforward to modify our formalism to ac commo- 
date other re asonable assumptions, as discussed in |Madavj 
et al. (1999) .) As the dust density declines with distance 
as R , the light "echo" observed at time < b = t(l + z) 
will be scattered by dust concentrated in a ring located at 
the intersection of the sphere R — i? su b and the paraboloid 



R x 

c 



(2-2) 



where fi = cos 9 (see Fig. [I]). In our model, it is adequate 
to ignore a finite re-processing time and the radial distri- 
bution of the dust. The dust only has to survive for a time 



From Eq. ( |2-3D it is clear that the only time dependence 
comes from the angular dependence of the escape prob- 
ability, dV sc /d£l; F ac (tob) is simply a step function for 
isotropic s cattering. W e adopt a Henyey-Greenstein func- 
tion (e.g. White 1979| ) to describe the differential cross 



section for dust scattering: 



da 



(A' 



dn 



[1 + W - 2</^] 3 / 2 



(2-4) 



with (fi) = 0.5 ( [Draine fc Lee 1984j )- We then use Eq. 
( |2-4| ) to compute dV sc ( A, JlJJdJl numerically for a slab-like 
dust cloud. The results for different observer angles (w.r.t. 
to the slab normal vector) and two different values of the 
total extinction, r = r a b s + t sc (measured at A = 0.3 (im) 
are shown in Fig. ||. The differential escape probabil- 
ity is normalized so that the integral J 4jr ^-(A,/i)rfO is 
equal to the escape probability from the dust cloud. Fig. 
|2| shows that at low optical depth r .3 = r(0.3 /xm) 3, 
the echo should have a similar color to the OT whereas at 



3 



larger optical depth, the echo will be much redder due to 
absorption. 

To illustrate that our calculation of the escape probabil- 
ity is not overly simplified (though it ignores wav elength 
dependence of the functional form for da/dfl, (White 
19790 in Fig. @ we show one curve (thin long-dashed line) 
computed using da/tKl oc 1 + 2/i + /i 2 . This expression 
gives the same v alue of (/i) but is less strongly peaked at 
fj, = 1 than Eq. (f^J). The resulting dV sc (X, /z)/dfi is very 
similar to what we use in our calculations. 

The angular dependence of the escape probability is 
exhibited in Fig. for a dust cloud with To. 3 = 7 and 
three values of the incident photon wavelength. Note that 
<fP sc (fi)/dft remains relatively flat for <, 6 SC ~ 20° and 
decreases exponentially at larger angles. 
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Fig. 3. — The differential escape probability is plotted as a 
function of for a dust cloud with the optical depth for extinc- 
tion to. 3 = 7. The results are shown for three different values of the 
incident photon wavelength. 

2.4. Infrared Echo 

Hot dust will also emit an isotropic infrared echo due to 
thermal emission from dust at the rapid sublimation tem- 
perature ~ 2300 K, peaking at an observed wavelength 



A ~ 2 (1 + z) /im. Waxman fc Draine (1999) argue that 



only the photons in the 1 — 7.5 eV range will contribute 
to dust heating. For To. 3 *~ 7 the absorption efficiency for 
photons in this energy range is > 0.8; and moreover, such 
photons are likely to carry a considerable fraction of the 
total OT emission. Therefore the integrated infrared flux 
is 

E _ LAt OT c 1 
FlR ~^DlRjlT^- (2_5) 

3. COMPARISON WITH OBSERVATIONS 
3.1. OT-Echo-Redshift Relations 
Ad opti ng our simple model of dust scattering, Eqs. 



( 2-1 , 2-2 ) allow us to relate the sublimation radius and 
OT power^ 10 47 L47 erg s _1 , to the observed echo delay, 
s. 



ob6 

^sub 
L47 



0.2C^C^ hfi (l + zy 1 pc, 
0.03(1 + z)- 2 Cr 2 C 2 - 2 (^ b . 6 ) 2 , , 



(3-1) 
(3-2) 



where C\ = (1 — /i)/0.06 allows for beaming or characteris- 
tic scattering angles different from 20°, and C2 = R/R S ub 
should be used if the dust is located beyond i? su b- 



For simplicity, we now suppose that the spectral index 
of the OT is a ~ 1. This is quite close to the spectr al in dex 
of the observed afterglows. We can then use Eqs( [2-3| ) to 
relate the R-band (0.65/im) echo flux density to the escape 
probability 



[0.65/im] - 0.4 



dn 



ob.6 ob,3 

CiC* 



(3-3) 



where the observed duration of the optical transient is 
10 3 A^ b T 3 s. Note the strong dependence on redshift which 
implies that accurate measurements of both the optical 
transient and the echo flux could lead to a fairly precise 
redshift prediction. 

The ratios of the optical transient flux density, F®J^ = 
iv/ns(47r) -1 Z}^ 2 (l + z) -3 , and infrared echo flux density 
to the optical echo flux density are likewise given by 



F° T [0.65/tm] 
FE[0.65/tm] " 

F^ E [2(l + z)/im] 
F„ E [0.65/im] 



3000 



0.5 



( tE 



ob.6 



dP si 



(l + z), 



(3-4) 
(3-5) 



where / ns is the fraction of incident OT photons, emerging 
unscattered from the dust cloud. 



3.2. GRB 980326 



For GRB 980326, 



an excess R flux [0.65/im] 



loom et al 



0.4 /i jy was measured a time t\ 

1999). If we make the simplest assumptions, a_i ~ 
Qabs ~ C\ ~ C 2 ~ 1, then R suh ~ 0.3(1 + z) _1 pc and 
L ~ 9 x 10 45 (1 + z)~ 2 ergs -1 . Comparing the reported 
spectral slope, a ~ 2.8 of the putative echo to that of the 
afterglow (a ~ 0.8), we estimate that t .3 ~ 7 (c/ Fig. 
2). This, in turn, implies that dP sc /dD, in the observed R 
band - 0.2(1 + z)~ l and / ns ~ 0.05(1 + z) -4 (see Fig 



We can then use Eq.(3-2) to deduce that Af^ r 3 ~ 3(1 +z 



and [o.65/im] ~ 200(1 + z) -10 fjjy. If z ~ 0.4, then 
the energy associated with the first optical measurement 
of the afterglow (F°t [0.65/im] - 10 /iJy after 0.5 d), suf- 
fices to account for the observed excess after 20 d as a dust 
echo. If z > 0.4, then the optical transient would have had 
to be present and create a larger fluence at earlier times. 
This is not unreasonable as the OT flux was measured 
to satisfy F OT oc <~ 2 . In view of the large number of 
simplifying assumptions that we have made, this estimate 
can only be regarded as illustrative. However it suffices to 
demonstrate that dust scattering is consistent with all of 
the available data. 

3.3. GRB 970228 

A somewhat similar story can be told for GRB 970228, 
whe re th e redshift, z = 0.695, is known (Djorgovski ct 
1999p . The earliest R-band measurement is ~ 30/iJy 
0.7 d after the GRB; and after ~ 30 d there red excess flux 
0.3/iJy was observed, with the spectral slope (a ~ 3.0) 



very similar to that seen in GRB 980326 (Galama et al 



4 



1999). For this object, again, within the uncertainties, the 
flucncc measured in the first stages of the optical transient 
is sufficient to account for the energy in the optical excess. 

3.4. Dust Origin 

In both examples above, the mass of dust required to 
produce an optical depth To. 3 ~ 7 with our simplest as- 
sumptions and assuming that it is spherically symmetri- 
cally distributed with respect to the GRB is ~ 0.1 M Q . 
This amount of dust could form in an expanding high- 
metallicity wind associated with an earlier stage in the 
evolution of the GRB progenitor as we have assumed in 
our simple model. Alternatively the dust might be associ- 
ated with a molecular cloud if GRBs are associated with 
massive star formation or a molecular torus should they 
be located in obscured galactic nuclei. 

4. DISCUSSION 

In this letter we present an alternative explanation for 
the reddened excess emission observed in GRB 970228 
and GRB 980326, which we attribute to dust scatter- 
ing of the early-time, afterglow emission. This sce- 
nario is predictive enough to be confirmed or ruled 
out with observations of future GRBs. In particu- 



lar, in contr a st to the supernova e xplanation ( Bloo 



Reichart 1999| J7 



1 



et al. 1999| ; palama et al. 1999| ; 

tEe" excess emission is due to dust scattering, then its 
properties will depend o n the luminosity of the opti- 
cal transient. HETE II flhttp: / /space.mit.edu/HETE/j ) 
scheduled to be launched in early 2000 and Swift 



(http:/ /swift, gsfc.nasa.gov/homepagc. html), scheduled 
for 2003, should provide real-time localization of GRB 
X-ray afterglows with sufficient precision to permit faster 
follow-up and better measurements of its total fluence. 
Infrared observations may discover th e expected thermal 
emiss ion from hot subliming dust (c/ Waxman & Drainc 
1999). In fact dust emission might be the correct explana- 



tion for the "near-IR" b ump seen in the s pectrum of the 
GRB 991216 afterglow (|Frail et al. 2000|). Note that as 



most GRBs are at redshifts <; 0.5, 3 /im, (as as opposed 
to the more common 2 /im) photometry may be necessary 
to see this emission. 

In those GRBs, where it is also possible to measure a 
redshift, the the simplest model of dust-scattering is over- 
contrained and therefore refutable. Beaming and dust 
inhomogeneity introduce additional uncertainty but such 
models may also be excludable. For example, if ROTSE 
dAkerlofet al. 1999| ) were to detect another optical flash in 
a GRB as luminous as that seen in GRB 990123, which had 
an isotropic luminosity L ~ 10 51 erg s _1 , then dust should 
be physically sublimed out to a distance i? su b ~ 100 pc 
along the line of sight. Unreasonably large beaming would 
then be required to explain a dust echo with a delay of 
only a few weeks. Alternatively, if the radio light curve 
in an afterglow tracked the optical light curve, then this 
would be incompatible with both dust scattering and a 
supernova. 

A further prediction of the dust echo model is that, un- 
less the dust and OT are both arranged axisymmetrically 
with respect to the line of sight, we expect there to be lin- 
ear polarization associated with dust echos and this may 
be measurable in bright examples. (1.7 percent polariza- 
tion has been reporte d in the optical trans ient associated 
with GRB 990510 by |Covino et al. (1999)| but this is un- 
likely to be due to scattering.) 

In conclusion, we have demonstrated that dust scatter- 
ing can account for the excess optical emission observed in 
the afterglows of two GRBs as an alternative to an under- 
lying supernova explosion. Future observations should be 
able to rule out or confirm this explanation. 
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